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1. Introduction 

In several fields of nuclear and condensed matter physics, as well as in astrophysics 
and plasma physics, there exists experimental evidences which suggest strongly the 
necessity of introducing new statistics, different from the standard Boltzmann-Gibbs, 
Bose-Einstein and Fermi-Dirac ones. 

Typically, generalized statistical distributions arise in systems exhibiting long-time 
microscopic memory, long-range interaction or fractal space-time constraints PIEIISI- 
Solar neutrinos |3], high energy nuclear collisions [S], cosmic microwave background 
radiation [HI IZI are some typical phenomena where generalized statistical theories has 
been applied. Several models of g-deformed statistics have been used in the study 
of Bose gas condensation [SI El, phonon spectrum in ^He jTU] as well as in the study 
of asymmetric XX Z Heisenberg chain ^1]. Moreover, quantum g-statistics has been 
introduced in the blackbody theory [T21 CH] where, based on an equilibrium statistical 
approach, a deformed Planck's law has been derived. 

Recently, in jHj, a generalized kinetic theory has been proposed by Rossani and 
Kaniadakis, with the purpose to deal, at a kinetic level, particles obeying to a generalized 
statistics. Following this idea, generalized kinetic theories of electrons and photons ^H] 
as well as electrons and phonons (THj have been proposed. 

In this paper we present, as a natural continuation of the previous works, a kinetic 
theory of interacting atoms and photons obeying to a very general statistics. 
We recall that in the last years a classical kinetic approach in the study of this 
dynamical problem has been proposed for the case of two energy levels |17| ITS] and 
further generalized for multi- levels atoms and multi- frequencies photons 1201 • The 
most remarkable feature is that this approach allows a self-consistent derivation, at 
equilibrium, of Planck's law. 

The physical system we deal with is constituted by atoms {i = 1, 2, ■ ■ ■ , N) with 
mass m, endowed with a finite number n of internal energy levels = E-^ < < 
■ ■ ■ < E^, with transitions from one state to another made possible by either scattering 
between particles, or by their interaction with a self-consistent radiation field made up 
by photons p.. with intensity L., i < j = 1, 2, ■ ■ ■ , n, a.t the n{n — l)/2 frequencies 
uj.. = E^ — E. (we adopt natural units h = c = 1). 

According to the relevant literature [21], we assume that the following interactions take 
place: 

a) elastic and inelastic interactions between particles, 

b) emission and absorption of photons. 

We restrict ourselves to the most common and physical interaction mechanism between 
particles: 



where i < j. The most general reaction ^Hj does not introduce significant improvement 
with respect to the present simplification. 
Gas-radiation interaction processes include: 



A, + A- ^ A, + A, 
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(1.2) 



(1.3) 



(1.4) 



Following Ref.s ^3 Cnj, for these emission and absorption processes we assume that: 
the outcoming atom has the same velocity as the incoming atom (that is we neglect pho- 
ton momentum with respect to atom momentum); photons are spontaneously emitted 
isotropically; photons which derive from stimulated emission have the same directions 
as the incoming photons. 

Finally, gas-radiation interactions are modeled by means of the Einstein coefficients |22j : 
for both absorption and stimulated emission, a^^ for spontaneous emission. 

The paper is so organized. In the next Section 2, we introduce the kinetic equation 
for the physical system of interacting atoms and photons. In Sections 3 and 4 theoretical 
results on equilibrium of the deformed model and its stability, by means of Lyapounov's 
theory, are given. The mathematical results are discussed on a physical ground, and 
connections with thermodynamics are pointed out. Two particular situations are studied 
in Section 5 and 6: photons in a background of atoms and atoms in a background 
radiation. In Section 7, we drop the statistics generalization for atoms but keep the 
generalization for photons. The results of two numerical simulations are presented. In 
the first simulation, the system initially in equilibrium, is disturbed, starting from time 
t = 0, by a fiash of monochromatic photons injected by an isotropic external source. In 
the second simulation, an isotropic source of monochromatic photons, constant in the 
time, is inserted starting from time t = 0. Section 8 deals with the conclusions. Finally, 
in the Appendix A we present an overall sum up on the kinetic equations for a system 
of atoms and photons obeying to the standard statistics while, in the Appendix B we 
recall briefiy two generalized statistical distributions used in the numerical simulations. 

2. Generalized Boltzmann equations 

In this section we introduce the kinetic equation describing a system of atoms and 
photons obeying to a very general statistics. For sake of clarity we begin by introducing 
the Uhling-Uhlenbeck equation (UUe) in order to fix the notations used in the following. 
The UUe is a quantum extension of the Boltzmann equation Let us consider a 

system of particles with mass m which obeies the Bose-Einstein or Fermi-Dirac statistics. 
These particles interact elastically by means of binary encounters. In this situation the 
UUe reads as 




gaig,C){fiv')fiw') l + A/(^) l + A/(^) 



> . (2.4) 
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- f{v) f{w) [l + A f{v')] [l + A f{w')] I dw dn' , (2.1) 

where A = 1 for bosons and A = — 1 for fermions. In Eq. ()2.H) f{v) = /(a?, t) is the 
dimensionless distribution function, a{g, Q is the cross section, with g = \v — w\ the 
relative speed. The post-colhsion velocities are given by 

v' = ]^{v + w + gn') , w' = ]^{v + w~gQ,') , (2.2) 

where v and w are the velocities of the incoming particles and cos^ = ^ ■ f^', with 
Q, = {v — w)/g the unit vector in the direction of the relative speed. Finally the two- 
dimensional unit sphere S'^ is the domain of integration for the unit vector fi'. 
In view of applications to particles which obey a more general statistics, we introduce 
the following substitution in the collision integral 

f{v) ^ if[f{v)] , 1 + A f{v) ^ ^[f{v)] , (2.3) 

where (p{x) and ip^x) are non-negative functions supposed to be continuous as their 
derivatives. The limit values for a; — > are subject to precise physical requirement. 
For instance, V'(O) = means that no transition occurs when one of the initial states is 
empty, ipi^) = 1 means that no inhibition occurs when one of the final states is empty. 
According to [12], we assume 

_d_ f i^{x) 

dx \ip{x) 

This is trivially true for bosons and fermions. In general this assumption is justified a 
posteriori since it assures uniqueness and stability of equilibrium. 
The generalized UUe now reads as follows 

1^ + ^ • V/ = j^^^^j^g. C)|^[/(^')] ^[/M] ^[/(^)] ^[/(^)] 

- ^[f{v)]v[f{w)]^lj[f{v')]^l,[f{w')]^dwdn' . (2.5) 

In the following we specify Eq. (j2.5|) for a system of interacting atoms and photons. 
Let us introduce the characteristic departure and arrival functions jEj: v?(/J and '?/'(/J 
for atoms, as well as ^{P^) and '^{P^) for photons, where atoms A^, endowed with energy 
level El,, are described by means of the distribution function /^(a;, t), and photons 
p.., at frequency u.., are described by means of the radiation intensity E^ix, f2, t). 
The distribution function of atoms Ai obeys to the following system of generalized 
Boltzmann equations (see Appendix A) JH] 

yi + v.Vf, = J,{v) + JXv) , (2.6) 

where in the right hand side of Eq. ()2.6p J^^v) describes the contribution due to the 
atom-atom interaction and Je{v) the contribution due to the atom-photon interaction. 
Explicitly, the collision integral, which takes into account the inelastic and elastic atom- 
atom collision, is given by 

^ (-) = E E + E E ^ (-) + E E + E E ^ (-) , (2.7) 

j i<j j i<j k j>£ k i<£ 
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J>) = j^^^s^sCig, C)|^[/.K)]^[/.K)]^[/.W]^[AM] 

- ^[f^{v)]v[f^{w)]ij[fXvimf,{wl)]^dwdn' , (2.11) 

where cx^^ and a"' are the cross sections for forward and backward reactions describing 
elastic and inelastic interactions. They satisfy the following microreversibility 
relationships [12] 

g'o';^{g, n') = {9Zfa;^{g-, n') , g^a^ig, U') = (4)V;;(4, n') . (2.12) 
The post-collision velocities are defined by 

vf^ = ^iv + w + gf^n') , njf^ = l(v + w-gf^n') , (2.13) 

with 



gf, = \ 9' ± - (E^ - EJ . (2.14) 
^ V m 



The four contributions to J),{v) correspond to the cases in which plays the role, in 
the reaction scheme (jl.ll) . of A^ in the r.h.s., A^ in the l.h.s., A. and A, respectively. 

Differently, the collision integral J^v), which accounts for the gas-radiation 
interactions, is given by 

where 

Uv, ^) = |^[/.(*^)] num mi-^)] - Mi^)] mAm mi-^)]] ■ (2.16) 

By taking into account all the energy levels higher than i, the loss term is due to absorp- 
tion, while the gain term is due to spontaneous and stimulated emission. The situation 
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is reversed when we consider all the energy levels lower than i. 



Finally, the kinetic equation for photons pij reads: 
dl 

+ =u..j..(n) , 



where 



dt 



j-.m 



J- .(v, fl) dv 



JR3 



(2.17) 



(2.18) 



Again, the gain term is due to spontaneous and stimulated emission, while the loss term 

is due to absorption. 

It is easy to see that, by posing 

^[/,(t;)]^l, (2.19) 



ij 



a. 



(2.20) 



Eq.s ()2.6p and ()2.17|) reduce to the standard kinetic equations for atoms and photons 

[111201. 

From now on, for sake of simplicity, we adopt the notations V'^('y) = ip[f^{v)], 

^,{v) ^ ^[f,{v)], ^ and ^^^{n) ^ w.M]- 

3. Equilibrium 

Notwithstanding the generalized equations are more complicate with respect to the clas- 
sical kinetic equations, many of the methods used in the standard kinetic theory are 
still applicable. This is in particular true for the study of equilibria and their stability. 



Lemma 1. For any given arbitrary smooth functions 7^(1^) and F-^ (f2) the following 
two relations hold 



q a (g, C) 



{v^) (uj+) V, {v) {w) - {v) Lp^ {w) (v+) (?i;+) 
7,(«) + 7,(^) - 7fK^) - 7.(^^,^)1 dvdwdQ' , 



(3.1) 



IR3 



3 i-<3 



52 



EE 

3 i-<3 



a. 



{v) {n) {v) - {v) {SI) i,^ [v] 



dv dfl 



(3.2) 



Proof. The proof for the generalized case follows the same steps which can be 
found in Ref. ^Hj for the standard theory. In the following we present the main points 
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recalling to the relevant references for the details jTOl |211 12S] • 

Eq. ()3.1|) follows from the definition of the collision integral J^iv) given in Eq. ()2.7p . 
The relevant addends may be grouped together, with suitable interchang of the indices 
i, j, k and £, to form the collision term appearing on the r.h.s of Eq. p.ip . In particular, 
for all the J^^^^{v) and the microreversibility relationships ()2.12j) are used as well 

as the relation 

dv dw dn' = dv^ dw^ dn , (3.3) 

where the Jacobian of the transformation in Eq. ()3.3p arises from the definitions ()2.13|) . 
In the same way, Eq. ()3.2|) is obtained from the collision integrals J^{v) and J^^i^) 
given in Eq.s ()2.15p and ()2.18|) . after suitable interchange of the indices i, j and i. • 

In the space homogeneous case, where all the dynamical quantities and J.^ are 
functions only of the time, equilibrium is defined by 

Let us introduce the following functional: 

After using the kinetic equations (|2.6|) and (j2.17j) . Eq. (|3.5|) becomes 



^ = E / In ^ J^{v)dv 



+ ^.W^^ + EEX.ln(y Wdn, (3.6) 

and by applying Lemma 1, with 7^ = ln(y9^/?/'^) and T^. = ln($.^/\I'^^.), we obtain 



X In — ^ dv dw dn' 



+ E E % / , , b {v)^^^{n)ijXv)- V. (^) m 



X In ^ ,cl^ , ' 3.7 

which is a quantity manifestly not positive since (A — B) XwiBj A) < 0. 

Proposition 1. Condition (j3.4p is equivalent to the following couple of equations 

v^{v)^^^{n)^X^) = ^X^)^jn)^^{v) . (3.9) 
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Proof. Taking into account the definitions of the coUision integral given in Eq.s ()2.7|) . 

and (jTTR|l . we observe first that Eq.s (jSSl) and imply J^v) = J^v) = 

J-^(r2) = and thus, from the kinetic equations ()2.6|) and ()2.17p . we obtain condition 
(III. 

On the other hand, from Eq. ()3.5|) by using Eq. ()3.4j) . it follows that V = 0. But, since 
in Eq. ()3.7|) both the integrands are never positive, condition V = implies Eq.s ()3.8p 
and dSSl). • 

In the case i = j, Eq. ()3.8j) gives 



where T is the absolute temperature of the whole system of atoms and photons. (Here 
and in the following we pose k-Q = 1). 

For the case i j, from Eq. ()3.8p and by using Eq. ()3.1H) . we get, as additional 
condition, that all the chemical potentials are equal: fi- = fi. = fi. 
Differently, from Eq. ()3.9|1 . using Eq. ()3.11|1 . follows 



which is the generalized version of Planck's law. In Appendix B we write explicitly Eq. 
(|3.12p for two particular statistics. 

Let us to observe that, due to the monotonicity of both (pjipf> and ^../"if.., Eq.s 
()3.11|) and ()3.12|) give an unique equilibrium solutions for f^{v) and L.{fl), respectively. 

4. Stability 

In order to study the stability of such equilibrium solution let us introduce the following 
functional 




(3.10) 




(3.11) 




(3.12) 



L = Ha + H, 



p ' 



(4.1) 



where 




J i'<3 ^3 



(4.3) 



(4.2) 
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with 



df, 
dl 



- In 



(4.4) 
(4.5) 



Remark that, since ^Jip^ and have been assumed to be monotonically increas- 

ing, Ha (/J and Hp{L-) are convex functions of their arguments. 

From P3] we know that 5* = —L is nothing else but entropy density for the present 
physical system and it is the sum of two contributions: gas entropy and photon entropy 

EH El ESI. 



E 



Lemma 2. The condition 

'dHAif, 



(4.6) 



where * means at equilibrium, is equivalent to energy conservation E{f^, L.] 



Proof. Energy can be written as 
E = Ea + Ep , 



where 



Ea = / SaM dv 



E„ 



S2 



with 



Sp{n) dn 



1 



^p = EE4(f^)- 

By means of Euler's theorem we can write 
'OSa ■ 



dn 



3 «<j 



Differently, by using Eq.s (pTTTI and (pTT^ . from Eq.s (|01) and (031), the 
relationships are easily obtained 



dE 



1 

2^* 



^ dSp 
T* [dl 



(4.7) 

(4.8) 
(4.9) 

(4.10) 
(4.11) 

(4.12) 
following 

(4.13) 
(4.14) 
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From Eq. (j4.12j) . by taking into account Eq.s (j4.13p - (j4.14j) and particle conservation 
J2e J fe dv = J2e I f* dv, one easily realizes that E = E* implies Eq. fl4.6|) . • 



The main result with respect to stability can be summarized as follows: 



Proposition 2. The functional L is a strict Lyapounov functional for the present 
dynamical system: the unique equilibrium is asymptotically stable, and any initial state 
will relax to it asymptotically in time. 



Proof. First one easily proves that 
dL 



dt 



V<0 



(4.15) 



with dL/dt = only at the unique equilibrium position. In fact, from the definition 
dUD) it follows 



dL , , 

~dt^\ 



d-HAif,) df, 



d^ + EE- 



1 f dn,{L,)dL, 



dt 



dfl 



(4.16) 



V df, dt . 

and taking into account Eq.s ()4.4j) and ()4.5|1 . from the definition follows Eq. ()4.15|1 . 
On the other hand, by taking into account Lemma 2, one has 

r r * TT TT* I U TT* 

L — L — Ha — -Ha + -"p ~ -"p 

= E /3^A(/jrf«-E/3^A(/;)rf^^ 



dn 



3 *<J *J 3 i<3 »J 



df, 
dfl 



(/. - /;) dv 



3 i<j ^^i -^^^ 



EE- 



j i<3 ^^J •'^^ 



9/, 



(L.-r)dn 



1 



E ^a(/j rf« + e E — L ^p(4) 



3 i<3 ^'J 



where 



nA{f,) = nAif, 
7?p(/J = 7^p(4, 



^a(/;) + 



. df, 

' dn,{L;, 

dL 



a - /;) 



(/ -F' 



(4.17) 

(4.18) 
(4.19) 



Since both LCa and Tip are convex, we have LIa > and Tip > 0, and then we can 
conclude that L > L* and attains its minimum, L = L*, only for the unique equilibrium 
position. • 
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In the next two Sections we inquire on equilibrium and stability in two particular 
cases: photons in a background gas and atoms in a background radiation. 



5. Photons in a background gas 



An usual assumption in radiation gasdynamics is that the relaxation time due to atom- 
atom interactions is much quicker than the one due to gas-radiation processes. Thus, in 
the kinetic equation for photons we can fix the distribution of atoms as an equilibrium 
function /* at a certain temperature T. 

In order to study equilibrium and its stability, we introduce the following functional 



. . ■ U J '?2 



In 



exp 



T 



'dl 
dt 



(5.1) 



where o means that has been substituted by /*. 
After using the kinetic equation ()2.17|) Eq. (j5.H) becomes 



In 



52 



# """P (v# 



■I An] 



dfl 



(5.2) 



and using Lemma 1 with all the 7^(1^) = and T.^. (f2) = ln[($.^./\l/.^.) exp(a;.^./T)] we 
obtain 



Cp = EE 



a. 



52 



In 



exp 



T 



Finally, observing that 

and taking into account Eq. ()3.1H) . the quantity Cp becomes 



dvdfl . 

(5.3) 

(5.4) 



Cr. 



EE 



a. 



52 



X m — , , , — , , , di 2 



(5.5) 



which is manifestly not positive. 

From usual arguments (see proof of Proposition 1), the equilibrium condition d Ljdt = 
is equivalent to 

cp*iv)<ij^^{n)^ij*{v) = ^*{v)%in)r^iv) . (5.6) 

Again, after using Eq. (j3.11|l . from Eq. ()5.(j|l we find the generalized Planck's law ()3.12j] 
for photons. 

In order to investigate the stability of such equilibrium, let's introduce the following 
functional: 

Er. 



(5.7) 
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where Hp is defined in Eq. (j4.3|) and Ep is given in Eq. (|4.9p . 



Proposition 3. Lp is a Lyapounov functional for the present problem. 



Proof. First of all we find 
dLr, 

Cp<0, 



dt 



thus Lp is a decreasing function. In fact, from the definition ()5.7|) we have 

'dHpiL 



d Lp ^ x ^ -k ^ 1 



UJ.. 



dL 



T 



'dLV 



(5.9) 



and after using Eq. ()4.5j) we obtain Eq. ()5.8|1 . 
Moreover, since 



ij_ 

' T 



V dL / 

as it follows from Eq.s (|3.12p and (j4.5j) . we can write 
Lp - L*p = Hp - Hp + — (^Ep - Ep) 



(5.10) 



/* dn 



1 



i i<j ^^J •'^^ 



EE 



J i<j ^ij ■'^^ 



np{L^)-np{r) 

7?p(/Jdf2, 



dL. 



L -I* 

13 13 



dfl 



where Hp is given in Eq. ()4.19|1 . Because Hp is convex we conclude that Lp 
is Lp attains a minimum at equilibrium. • 



(5.11) 



> L*p, that 



Taking into account of the definition of Lp it is easy to realize that Eq. ()5.8|) is 
equivalent to Clausius inequality 



dSr, 



1 dEr, 



where Sp = —Hp. From Lp — L* > follows that the quantity Sp = Sp — Ep/T is always 
increasing and attains a maximum at equilibrium. 



6. Atoms in a background radiation 

Suppose now that temperature is high enough so that we can consider photons as an 
equilibrium background. This means that, in the kinetic equations for atoms, we can 
fix the distribution of photons as an equilibrium function /* at a certain temperature 

r. 
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In order to study equilibrium and its stability, we introduce the following functional 

1/1 o „ Ml (dfy 



dt 



dv 



(6.1) 



where jj means that L. has been substituted by I* . After using the kinetic equation 
(EH), Eq. dniH) becomes 

1 n 



E 



In 



iR3 



exp 



T V2 



mv + E^, — 



dv . 



(6.2) 



Using Lemma 1 with 'y^{v) = ln{{ipjilj^) exp [{mv'^/2 + E^ — fi) /T]} and all 
0, Eq. ()6.2|) becomes 



= yyy g^ig.o 



X In ■ 



X exp 



-E. 

3 l_ 

T 



m 
2T 



X In 



[^^,(^) r,(f7) ^^(i;) - cfXv) r,(17) V^^.(t;)] 



exp 



dv dn 



_^^.(^)^,(t;) V T 
and taking into account the energy conservation and Eq. (|3.12|) 



ui. 



E. -E. = -m 
J . 2 



V. 



( +\2 2 2 



(6.3) 



(6.4) 



we obtain 



g cr (g, C) 



= EEE 



X In 



<^£(^+)<^.("^+)V',(^)^,(^) 



+ EE 



X In 



3 «<J 

¥.,(«)$* (17) V^,(t;) 



(iv dft 



(6.5) 



which is manifestly non positive. 

Following the same steps given in Proposition 1, we obtain that the equilibrium condition 
dfjdt = is equivalent to 
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The first Eq. gives 

ln^ = -^ i^-mv' + E^-li^ , (6.8) 

while the second Eq. flfj.Tj) . after using Eq. (jfi.Sj) gives T = T. 

In order to study the stabihty of this equihbrium we introduce the following 
functional 

Ea 



La = Ha 



T 



(6.9) 



with Ha and Ea given in Eq.s (|4.2p and (j4.8|) . respectively. 



Proposition 4. La is a Lyapounov functional for the present problem. 



Proof. First we find 
dLA 



dt 



Ca < 



so that La is a decreasing function. In fact, from the definition ()6.9|) we obtain 



dLA 
dt 



E 



+ — - m + 



dt 



dv 



(6.10) 



(6.11) 



df, ■ T \2 

After using Eq. ()4.4j) and taking into account for particle conservation we obtain the 
relation (jOTl|l . 
Moreover, since 

fdUAif, 



df, 



(6.12) 



as it follows from Eq.s (|3.11|) and (j4.4p . by accounting for particle and total energy 
conservation we can write 

La - L\ = Ha - HI + — {Ea - E^) 

=i:L«A(/,)*-i:4«A(/;)*+^i:t (i"-^+^,-.) (/,-/;) 



iR3 



^a(/j-Ha(/; 



df, 



dv 



(6.13) 



with Ha given in Eq. (j4.18|) . Due to the convexity of Ha we can conclude that La > L*p^ 
and La attains a minimum at equilibrium. • 

Finally we can verify that Eq. ()6.10|) is equivalent to Clausius inequality 

where Sa = —Ha- From La — L^ > follows that the quantity iSa = "Sa — Ea/T is 
always increasing and attains a maximum at equilibrium. 
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7. Numerical Simulations 



For the purpose of numerical calculations, we consider a homogeneous and isotropic case 
and keep the generalization for photons only. Atoms are treated as classical particles. 
This means that substitution ()2.19|1 is performed, while ()2.20|1 is not applied. A couple 
of different cases are considered for the functions $ and \1'. 
By integrating Eq. ()2.6p over dv we obtain 

^ = E(^« +4vri?J - E(^« +4^^J , (7.1) 

where 

N^= f J,{v)dv, (7.2) 

is the number density of atoms at level i and 

S.,=T.L , ga2{g,0[fXv)f^{w)-fSv^)f\{w^)]dvdwdn', (7.3) 

A, =«.,[Ar^^(Jj-Ar^$(/J], (7.4) 

are the source terms due to the elastic/inelastic and emission/absorption interactions, 
respectively. 

By multiplying Eq. ()2.6p by m v^/2, by summing over ^ and by integrating over dv we 
obtain 

^N^ = yy^u..s.. , (7.5) 

with the total number density of the mixture N = J^e^^e a constant, and temperature 
given by 

^ TtF , (7.6) 



3N 
where 



]P = mY,f (v ® v) f^{v) dv , (7.7) 
e Jib? 



is the stress tensor. 

Differently, by integrating Eq. ()2.17j) over dQ. follows 
dl 



dt 



uj B . . il.i 

11 !J V 



Suppose now that elastic atom-atom interactions prevail over the inelastic and 
gas-radiation ones (consistently with the usual local thermodynamic equilibrium 
approximation). Thus, according to the zero-order Chapmann-Enskog approximation 
we can calculate 5*. . by means of the Maxwellians 

It is found 

(7.10) 



S.. = G 

'J I] 
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with 



G,=j:N,\iT), (7.11) 

e 

%(T) = (^Y [ ^ ^ x) exp {v' + w')] dvdwdn' . (7.12) 

Remark that the equihbrium condition is equivalent to 5 = B =0, from which follows 
the relation 

K exp (^) = iV, exp ( ^ ) . (7.13) 



T 

Taking into account the total number conservation, we obtain the atoms distribution 
function 

N, = -exp(--^] , (7.14) 



with 



T 



Eexp(-|f) . (7.15) 



We observe that Eq.s ()7.14|) and ()7.15|) assume the standard expression for the 
equilibrium distributions of the atomic levels. Notwithstanding, because the generalized 
statistics is kept for the photons, the equilibrium temperature is different respect to the 
case in which also photons obey to the standard Bose-Einstein statistics. Consequently, 
the values of the population levels A^^ are numerically different. 

Eq.s fl7.1|) . ()7.5p and ()7.8|) constitute a close system for the unknown quantities N^, 
T and L.. They can be integrated by means of numerical calculation, with suitable 
initial conditions. 



In the following, we consider a system of atoms with n = 2 levels and 
monochromatic photons with frequency cu^j. The system is initially in the equilibrium 
configuration at temperature T*. From condition B^^^ = follows the distribution I*^. 

Two different simulation are take into account. In the first simulation, a flash of 
photons is injected at the time t = 0, from an isotropic external source at temperature T'. 
The relaxation of the system to the new equilibrium configuration is depicted v.s. time 
in the following figures 1 and 2, concerning the results of two different generalizations 
adopted for the photons and described in the Appendix B. They take into account for 
an asymptotic inverse power law decay of the distribution function with respect to 

First we adopt the model proposed by Biiyiikkilig et al. ^71I2H1- The time evolution of 
the physical meaningful quantities A''^, N^, I^^ and T v.s. the deformed parameter q is 
depicted in figures (la), (lb), (Ic) and (Id), respectively. 

Then, we adopt the model proposed by Kaniadakis [7j. The time evolution of the 
physical meaningful quantities A^^, N^, I^^ ^^"^ ^ ^-S- the deformed parameter k is 
depicted in figure (2a), (2b), (2c) and (2d), respectively. 
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1a lb 




Figure 1. Plot (in arbitrary units) of time evolution of N-^ (la), (lb), (Ic) 
and T (Id) v.s. q, after injection, at i = 0, of a flash of photons obeying to the g-Bose 
statistics. 



In the second simulation, a constant and isotropic external photon source is inserted 
at the time t = 0. Again, both the g-deformation and the K-deformation are considered. 
The time evolution of the same physical quantities N^^, N,^, I^^ and T v.s. the deformed 
parameter q [figures (3a), (3b), (3c) and (3d)] and v.s. the deformed parameter k [fig- 
ures (4a), (4b), (4c) and (4d)] are plotted. A comparison between the g-deformed case 
and the K-deformed case leads to the following considerations. With the respect to the 
classical case, the more the system is K-deformed, the steeper the evolution is. On the 
contrary, the more the system is g-deformed the softer the evolution is. 
We remark that this opposite behavior between the g-deformed case and the K-deformed 
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Figure 2. Plot (in arbitrary units) of time evolution of N-^ (2a), (2b), /^^ (2c) and 
T (2d) v.s. K, after injection, at i = 0, of a flash of photons obeying to the K-Bose 
statistics. 



case is a consequence of the choice g > 1. Such a choice arises from the observation 
that, in the present problem, for values of g < 1 a cut-off in the energy spectrum is 
required, (see Eq. (jB.ljl in the Appendix B). 

Moreover, from the second simulation we observe that, as well-known for the non- 
deformed case, in presence of a constant pumping of photons, the two number densities 
reach asymptotically, as t ^ oo, to the same value. Such a feature is preserved also in 
the deformed cases but with a different time-scale. 
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3a 3b 




Figure 3. Plot (in arbitrary units) of time evolution of N-^ (3a), (3b), /^^ (3c) and 
T (3d) v.s. q, after the introduction, at i = 0, of a constant photons source obeying to 
the q-Bose statistics. 



8. Conclusions 

In order to improve tlie knowledge of the generalized statistical mechanics and 
generalized kinetic theory, we have investigated, in the present paper, a system of 
atoms and photons obeying to a very general statistics. The kinetic equations has been 
obtained, in the extended Boltzmann picture, through the introduction of characteristic 
departure and arrival functions for atoms and photons. In the space homogeneous case 
we have studied the equilibrium configuration. Such equilibrium is given by a modified 
Planck's law for the radiation and by a generalized distribution function for the atoms. 
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Figure 4. Plot (in arbitrary units) of time evolution of N-^ (4a), (4b), /^^ (4c) and 
T (4d) v.s. K, after the introduction, at i = 0, of a constant photons source obeying 
to the K-Bose statistics. 



By means of Lyapounov's theory we have studied on the stabihty of this unique 
equihbrium configuration which maximize the entropic functional of the system. 

In the second part of this paper we have considered a homogeneous and isotropic 
case keeping the generahzation for photons only. Atoms are treated as classical particles. 
Under suitable assumptions on the relaxation times of the various interaction processes, 
according to the zero-order Chapmann-Enskog approximation, we have obtained a close 
system of macroscopic equations for the unknown dynamical quantities: the number 
density of atoms at energy E^, the temperature of the system T and the intensity 
/ . at frequency tu^j. 
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We have shown the results of some numerical simulations for a system with n = 2 energy 
levels. In the first simulation the system, initially at equilibrium at the temperature T*, 
is perturbed by a flash of photons injected by a external source at a different temper- 
ature T'. The relaxation of the system to the new equilibrium state for the physically 
relevant quantities N^, N^, T and has been plotted. In the second simulation the sys- 
tem, initially at equilibrium at the temperature T*, is disturbed by a constant source of 
photons. The time evolution of the same physically relevant quantities has been plotted. 



Appendix A. 

For sake of completeness, we report in this Appendix on the classical kinetic equations 
for atoms and photons, as known in literature. 

The distribution function of atoms A£ obeys to the following system of Boltzmann 
equations 

dfe 



dt 



v-Vf, = J,{v) + J^{v) 



(Al) 



where in the right hand side of Eq. ()A.1|) J^{v) and J^(t^) describe the contribution due 
to the atom-atom interaction and the contribution due to the atom-photon interaction, 
respectively. 

In detail, for elastic/inelastic interactions between atoms we have 

^(1) / N j-(2) / ^ , v^v^ t(3)/__\ , V^V^ yW , 



^ (-) = E E C (-) + E E C (-) + E E (-) + E E ^ (-) 

j i<j k j>l k i<£ 



(A.2) 



with 





j 






» = 


[ 


f" 

i]l 


» = 


[ 


(3) 


» = 


j 


ike 


» = 


f 


tj 

! a 

ei 


and 


H 

a are 



ga Ug^C) f,{v^) f^w^) - f,{v) f^{w) 



dw dfl' 



ga;^{g, () f.iw') - f,{v) fXw)\ dwdn' , 

g al {g, () [f\ {v' ) f\ {w;^ ) - f,{v) f,{w)] dw dQ' 
ga^ig, () If^v^J f.iw^J - f,iv) f,iw)] dw dn' . 



(A.3) 
(A.4) 
(A.5) 
(A.6) 



elastic and inelastic interactions. 
In Eq.s ffO|l - (|0|l we have 



= -(v + W + n,') 



w^ = - (v + w — q^ Q''^ 



Ig^ ± - (E, - 
m 



Q = - (v — w) , 

9 



g = \v — w\ 



cos ( = Q ■ fl' 



(A.7) 
(A.8) 
(A.9) 
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and the two-dimensional unit sphere S"^ is the domain of integration for the unit vector 
f2'. The four contributions to J^('y) correspond to the cases in which plays the role, 
in the reaction scheme (jl.lll . of A^ in the r.h.s., A^ in the l.h.s., A. and A-, respectively. 
For gas-radiation interactions, one can write 

li^) = E L n)dn-Y: f U^, n) dn , (a.io) 

where 

Uv, n) = fXv) + A, - A. /.(^) i^m ■ (A.ii) 

By taking into account all the energy levels higher than i, the loss term is due to 
absorption, while the gain term is due to spontaneous and stimulated emission. The 
situation is reversed when we consider all the energy levels lower than £. 
The kinetic equation for photons Pij with intensity /.^. reads '■ 

+ n-vK = ujJ^^{n) , (A. 12) 



dt 
where 



X(f7)= l^^.l^{v,n)dv. (A. 13) 



The gain term is due to spontaneous and stimulated emission, while the loss term is due 
to absorption. 



Appendix B. 



We describe briefly two models used in the numerical simulations. They take into ac- 
count an inverse power law decay of the distribution function with respect to energy 

miaisj. 

1) The first model that we consider has been proposed by Biiyiikkilig et al. [2ZII2H1- 
Relevant applications to the blackbody problem can be found in |12[ I13j. Preliminarily, 
we define the g-deformed exponential and logarithm 



exp^^xj 



ln„ 



X 



l-q 



X] 



(B.l) 
(B.2) 



l-q 

where = x for x > and = for x < 0. These g-deformed functions have been 
introduced in thermo-statistics for the first time in Ref. |29j . 
The characteristic functions $ and ^ are given by ^H] 

<l>(x) = exp In^(x) , (B.3) 
\E'('x) = exp ln„ (x) - ln„ ( ^^^^ — ) , (B.4) 



1 + x 
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where 1 < q < 2 and for g = 1 the standard statistics is recovered. For bosons we 
observe that the condition \l/(+oo) = +oo does not hold. 
The deformed Bose-Einstein distribution function is given by 

^1 



(B.5) 



2) The second model that we consider has been proposed by Kaniadakis [Zj 
Preliminarily we define the K-deformed exponential and logarithm 

1/k 



expj^j [X) 



VT 



KX 



X — X 



2k 



Remarkably the ^-exponential satisfy the property exp^^j(x) exp^^j(— a;) 
characteristic functions $ and \E' are given by J3] 

$(a;) = exp 



exp 



X 



(B.6) 

(B.7) 
l.The 

(B.8) 
(B.9) 



where < < 1 and for k = the standard statistics is recovered. In the case of 
bosons, the condition for $(+00) = +00 is fulfilled as well. The deformed Bose-Einstein 
distribution function is given by 



a. 



13. 



exp 



-\ -1 



- 1 



(B.IO) 
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